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We study the dynamical dielectric function of a two-dimensional hole gas, exemplified on 
[001] GaAs and InAs quantum wells, within the Luttinger model extended to the two lowest 
subbands including bulk and structure inversion asymmetric terms. The plasmon dispersion shows 
a pronounced anisotropy for GaAs- and InAs-based systems. In GaAs this leads to a suppression 
of plasmons due to Landau damping in some orientations. Due to the large Rashba contribution 
in InAs, the lifetime of plasmons can be controlled by changing the electric field. This effect is 
potentially useful in plasmon field effect transistors as previously proposed for electron gases. 
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I. INTRODUCTION 

Spin-orbit coupling (SOC), as one of the most 
important consequences of Dirac's theory for particles 
in solid state physics, is widely believed to play an 
important role for establishing electronic devices such as 
spintronic transistors.^^ By varying the single-particle 
properties of materials with sufficiently large spin- 
orbit coupling, typical response quantities such as the 
conductivity or the dielectric function can be controlled 
in an efficient way!^' The latter is not only necessary in 
understanding screening of extrinsic charged impurities, 
which in turn is important for transport, but also because 
of the existence of collective charge excitations known as 
plasmons 

Many analytical and numerical studies have been made 
in the last yea rs regarding the dielectric properties of 
electrorpl^jED] anc j no i e gag S y S t eml |LlltII] or promising 

materials like graphenePSiH] As has been shown in recent 
works, large analytical progress could be made in a two- 
dimensional electron gas (2DEG) including the effect of 
an asymmetric confinement, the Rashba SOC, and the 
contribution due to bulk i nversio n asymmetry (BIA), the 
so called Dresselhaus SOCp^lifllor in graphene including 
several types of spin-orbit interactions (SOIs) Due 
to the more complicated nature of the valence bands, 
the calculation of the dielectric function in hole gas 
systems turns out to be a formidable task even without 
SOIs. While for the three-dimensional case the free 
polarizability has been obtained for arbitrary frequencies 
and wave vectors within the axial modelp^l the solution 
for the two-dimensional case is not known so far. 

In this work, we use Luttinger's four-band model 23 
extended to the two lowest subbands including the 
lowest order Dresselhaus and Rashba contribution and 
discuss the dielectric functions of GaAs- and InAs-based 
quantum wells. As we will show, the plasmon mode 
exhibits a pronounced anisotropy being much stronger 
in GaAs than in InAs. By applying an electric field 
perpendicular to the 2D hole gas in InAs, giving rise to 
a tunable spin-splitting of the bands due to the Rashba 
coupling, it is possible to modulate the damping rate of 
the plasmons. 



This paper is organized as follows. In Sec. |TT] 
we introduce Luttinger's model including SOIs of 
the Dresselhaus and Rashba type and summarize the 



formalism of random phase approximation!^ In Sec. Ill 
the energy spectrum and the density of states of a GaAs 
and InAs quantum well is discussed. In Sec. |IV| the 
anisotropy of the plasmon spectrum and the influence of a 
finite electric field on the lifetime of charge excitations is 
studied and possible applications are described. Finally, 
in Sec.lV] we give a brief summary and a short outlook. 



II. THE MODEL 

A. Hamiltonian 

We start with the basic part of Luttinger's 
HamiltoniarPi for III-V semiconductors in three 
dimensions (setting li = 1 throughout this work) , 
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where c.p. stands for cyclic permutation of variables, 
Ji are the usual spin-3/2 matrices, and mo is the bare 
electron mass. 

Often Hq is simplified by setting 72 = 73 = 7 = 
(272 + 373) /5. The advantage of this so-called axial 
(in two spatial dimensions) and spherical (in three 
dimensions) approximation is the relative simplicity of 
the energy spectrum. In three dimensions the subsequent 
twofold spin-degenerated bands read 
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where m^/i — mo/(7i =F 2 7 ) is the heavy and light hole 
mass and h = (k x , k y ) the wave vector in the x-y plane. 
While these bands resemble that of an ordinary electron 
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TABLE I: Band structure and spin-orbit coupling parameters of GaAs and InAs. Taken from Reff 
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gas, the eigenstates are clearly more complicated^ due 
to their non trivial spin structure. 

In addition to the Hamiltonian in Eq. we take into 
account SOIs of the Dresselhaus type up to the third 
order in k (which indeed is the lowest order to describe 
the effect due to BIA)J2H 
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and the RashbgP^ contribution, 
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derived perturbatively by using the Loewdin 
approximation.^ While the former is always present 
provided bulk inversion symmetry is broken, as is 
the case for zinc blende structures, the latter arises, 
e.g., due to a suitable external electric field breaking 
structure inversion symmetry. Notice that the signs of 
the Hamiltonians Hq, Hp, and Hr are inverted such 
that all energies are positive. 

We now add a spatial confinement in the [001] 
direction and express the Hamiltonian in the basis of the 
subband functions tp n (z) = 2 / d sm(nn z / d) (n 6 N), 
with d being the width of the quantum well. It is 
worth noticing that both the contribution due to the 
Dresselhaus SOI as well as the one due to the Rashba 
SOI in the two-dimensional hole gas (2DHG) depend on 
the growth direction of the crystal.^ However, we left 
out the terms in Eq. Q which result in a dependence 
on the crystallographic orientation of both the applied 
external field and the momentum k since they can be 
considered as much smaller than r |^ 8,u |2ll 
The new matrix contains blocks of four-dimensional 
Hamiltonians, each given by the projection 

(p n \Ho + Hd + HR\(p m )- For k z this leads to 
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Often the off-diagonal contributions (n ^ to) are 
neglected. In this case the mode n is a good quantum 
number and ip n (z) is the exact eigenfunction in the z 
direction. However, this approximation is only justifiable 
if 73 <C 71 which is not the case for GaAs and InAs; see 
Table |H 



Dielectric function 



The charge susce ptib ility of the non interacting system 
can be expressed as^' 
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where the E\ are the eigenenergies and \X\) the n,-th 
components of the eigenstates (i = 1,...,4). The above 
expression contains a summation over all band indices 
and modes n«. The form factoi^Il F™ 1 ™ 3 (g), defined as 
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arises due to the spatial confinement in the growth 
direction. For the two lowest subbands only six 
independent form factors remain due to symmetry, i.e., 
FH F 2 l FH FH, FH and F%. While the latter 
two vanish, the four finite contributions, which have been 
calculated numerically, are shown in Fig. [T] 

The dependence on temperature and chemical 
potential is contained in the Fermi distribution function 
f{E). In this work, we assume zero temperature such 
that f(E) = 6[fi — E], with fi being the chemical 
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FIG. 1: (Color online) Independent finite form factors F^ni 



as defined in Eq. (p3| 
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potential. The general relation xo(?,^ w ) = [xoiq, 1 ^)]*, 
valid for all response quantities, insures that we can 
restrict ourselves to positive frequencies ui. 

Having found the non interacting charge response, 
the dielectric function in random phase approximation 
(RPA), i.e., the system's response to a scalar probe 
potential, immediately follows frorrP^ 

e(q,u) = 1 - V(q)xa{q,uj). 

V(q) — e 2 /(2eoe r q) is the Fourier transform of the 
Coulomb potential in two dimensions, eo the vacuum 
permittivity, and e r the background dielectric constant. 
Note that as the Hamiltonian is anisotropic, the dielectric 
function will be a function of both q = \q\ and the polar 
angle given by tan <fi q = q y /q x . 
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III. DISPERSION AND DENSITY OF STATES 
OF A QUANTUM WELL 

Before we start with the discussion of the dielectric 
function, we have a closer look on the energy spectrum 
and the density of states (DOS) 

A k 

The latter is related to the imaginary part of 
the polarizability through the Dirac identity, 
Im{l/(x±i0)} = Tn5(x), applied to Eq. 

In what follows, we choose parameters according to 
a GaAs and InAs quantum well grown in the [001] 




FIG. 2: (Color online) Top: Energy spectrum of (a) GaAs 
and (b) InAs for two different orientations (f>k = 0° (solid 
black line) and cj>k = 45° (dashed red), with tan^ = 
k y /k x . The dashed green line indicates the Fermi energy Ef 
corresponding to a carrier density of n = 5 • 10 14 m~ 2 . The 
inset in (a) shows the energies E^.2 and En around the T 
point. Bottom: Density of states of (c) GaAs and (d) InAs. 
In all cases, the electric field is set to E z = 1.5 • 10 7 V/m. 



FIG. 3: (Color online) The solid lines show the plasmon 
dispersion of GaAs for an angle orientation of <f) q = 0° (black) 
and (j) q = 45° (red). The dashed lines display the boundaries 
of the single-particle continuum for <j> q — 0° (black) and 
(j> q = 45° (red). The electric field is set to E z = and the 
hole density to n = 5 • 10 14 m~ 2 (Ef = 0.5 ■ E ). 



direction^ with a thickness of d = 20 nm, an electric field 
of order 10 7 V/m, and a hole concentration of n = 5 TO 14 
m -2 |28] -j^g Dan d structure parameters can be found in 
Table |H We restrict ourselves to the two lowest subbands 
(/ = 1,2) and neglect higher ones, as well as electronic 
and split off bands. Our numerical inspection shows 
that the dominant contributions in the dielectric function 
arise from the intraband and interband transitions with 
final states in the ground state light and heavy hole and 
in the first excited heavy hole band, respectively. Already 
the inclusion of the I = 2 light hole states does not lead 
to significant changes in the plasmon spectrum. Hence 
the influence of energetically higher bands such as the 
I = 3,4, ... subbands or the electronic and split-off bands 
will be even smaller. 

For convenience, we furthermore introduce the energy 
scale Eo = tt 2 Ji/ \2m, d 2 ) and use the notation that all 
wave vectors with a bar symbol have to be understood as 
dimensionless quantities measured in units of 7r/d; e.g., 
q = qd/ir. 



A. GaAs 

In Fig. [2j» the six lowest valence bands of GaAs are 
shown for an electric field of E z — 1.5T0 7 V/m. The solid 
black (dashed red) curve shows an in-plane momentum 
angle orientation of 4>k — 0° (45°), where tan<^ = k y /k x . 

Neglecting the off-diagonal contributions, the heavy 
and light hole energies of the nth subband right at 
the r-point can be approximated by E^/iniO) — 
n 2 (1 =p 272/71) E a . The two lowest bands, ^fti(O) = 
0.38 ■ Eq, are well separated from the others. However, 
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the four next states, Eya and En, are very close to each 
other, though they do not touch which can be seen from 
the inset in Fig.|2ja). Thus, in GaAs it is not sufficient to 
restrict to the ground state energies (n = 1) but one also 
needs to take into account the first excited subbandf^ 

For large enough momenta, k > 0.5, the anisotropy of 
the heavy hole bands is clearly larger than that of the 
light holes. Additionally, the Dresselhaus contribution 
leads to a spin splitting of the bands. As the Rashba 
part in GaAs is virtually negligible for realistic fields™ 
only minor changes occur in the spectrum once an electric 
field is turned on. 

The DOS obtained from Eq. Q is shown in Fig.^c). 
First of all, due to quantization in the z direction, 
the hole DOS is zero for energies smaller than Ehi{0). 
Second, as the energy spectrum shows singular points at 
E v h ~ 1-25 • Eq, see Fig. [2ja), we observe two nearby 
Van Hove singularities around E v jj- Without SOIs, spin 
degeneracy is recovered and both singularities merge. 
Finally, for large enough energies, E > 2 • Eq, the DOS 
remains roughly constant. 



B. In As 

The numerically calculated energy spectrum of InAs, 
including a finite electric field of E z = 1.5 • 10 7 V/m, is 
shown in Fig.[2jb). Compared to GaAs, the anisotropy in 
the spectrum of InAs is much weaker. As the anisotropic 
terms are proportional to 72 — 73, the reason for this is 
the smaller relative difference of 72 and 73, roughly being 
10 percent in InAs while in GaAs it is about 30 percent; 
see Table [I] Unlike GaAs, the Rashba contribution now 
dominates over the Dresselhaus part. While for zero 
electric field the band spin-splitting is small, a noticeable 
spin splitting, as shown in Fig. [2jb) , is mainly caused by 
the Rashba part. Another important difference between 
GaAs and InAs is that the r-point energy of the first 
excited heavy hole state, i.e., £^2 = 0.74 • E , is much 
smaller than that of the ground state light hole band, 
En = 1-81 • Eq. Hence the n = 2 subband cannot be 
neglected. 

The DOS of InAs shows features similar to those of 
GaAs; see Fig. |2jd). That is, the hole DOS is zero 
for energies E < E^i = 0.17 • Eq. The step at E^i 
describes the sudden occupation of the two lowest bands. 
Furthermore, singular points in the derivative of the 
energy dispersion at E « 0.70 • E and E ps En give 
rise to characteristic spin-split Van Hove singularities. 



IV. COLLECTIVE CHARGE EXCITATIONS 



where w q is the plasmon energy. For small damping rate, 
Eq. ^ can be approximated by 



Re{e(q,u q )} = 0. 



(9) 



In principle, one can solve Eq. ^ for arbitrary 
wave vectors. However, these mathematical solutions 
correspond only at sufficiently low damping to a physical 
resonance. A practical way to clarify this is to study 
the energy loss function, Im {l/e(q, u) + i0)}, which is 
relevant in experiments. 



A. GaAs 

In a recent worli^l using a multiband k ■ p model it has 
been demonstrated that plasmons in GaAs-based three- 
dimensional hole gas structures behave quite differently 
compared to their electronic counterpart. In particular, 
in three dimensions the plasmon energy always lies in 
the single-particle continuum, even for small momenta. 
Contrary to electronic systems, plasmons are therefore 
always Landau damped. 

In Fig. [3j we show the numerically calculated plasmon 
dispersion of a GaAs quantum well obtained from Eq. ^ 
for two different polar angles <j> q — 0° (solid black line) 
and <t> q = 45° (solid red). As mentioned above, the hole 



density is set to n — 5 • 10 



14 



(E F = 0.5 • E ); i.e., 



the two lowest bands are occupied. The black and red 
dashed lines indicate the boundaries of the single-particle 
continuum (SPC) in the respective orientation in order 
to see when plasmons acquire a finite lifetime. Region I 
is thereby caused by intraband transitions between the 
lowest lying valence bands (E^i), while the interband 
continuum II is due to transitions with final bands Ei l2 
and En, respectively. 

As can be seen from Fig. [3j the lower part of the 
SPC (region I) and the plasmon spectrum turn out to 
be clearly anisotropic. Furthermore, the (j) q — 0° mode 
becomes damped at smaller wave vectors compared to the 
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We now continue with the discussion of collective 
charge excitations . P lasmons are defined as zeros of the 
dielectric function)^ 



e{q,u} q ) = 0, 
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FIG. 4: (Color online) Energy loss function of GaAs for two 
different angle orientations, <f) q = 0° (black) and <f) q = 45° 
(red), and momenta (a) q = 0.6 and (b) q = 0.9. The electric 
field is set to E z = and the charge carrier concentration to 
n = 5 • 10 14 m~ 2 (E F = 0.5 ■ E ). 
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FIG. 5: (Color online) Angle dependence of the plasmon 
energy (solid line) and imaginary part of the dielectric 
function (dashed line) for GaAs with fixed q — 0.8. The 
electric field is set to E z — and the charge carrier 
concentration to n — 5 ■ 10 14 m -2 (Ef = 0.5 • Eg). 
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g = 45° solution. To see this, the energy loss function 
Im{l/e(qf,w + i0)} is plotted in Fig. gj For g = 0.6, 
see Fig. |4ja), in both cases distinct peaks occur in the 
loss function indicating truly coherent modes. For larger 
wave vector, q = 0.9 in Fig. Qb), the peak is smeared 
out in the 4> q = 0° direction, i.e., the mode is damped, 
while for <j) q = 45° it remains <5-like. 

In Fig. [5j the angle dependence of the plasmon energy 
and imaginary part of the dielectric function are shown 
for fixed q = 0.8. The energy (solid line) is redshifted for 
larger angles, up to a minimal value at <p q = 45°. The 
relative difference between the 4> q = 0° and <fi q = 45° 
energies is roughly 30 percent and thus of order of the 
difference of the hole gas parameters 72 and 73. The 
second (dashed) line in the plot shows the imaginary part 
of the dielectric function being proportional to the quasi 
particle lifetime. For cf> q < 29° the damping rate is finite 
with a maximal value at around <j> q — 18°. For 30° < 
4>q ^ 45° the imaginary part vanishes and the peak in 
the loss function becomes sharp. 

It should also be emphasized that for a given 
momentum, Eq. ^ leads to an additional solution with 
smaller energy. However, this solution does not fulfill 
the original condition of Eq. (|8| as the imaginary part 
of the dielectric function is not negligible. This is 
exemplarily illustrated in Fig. [6]ja) where the real part 
of the dielectric function (solid black) and the energy 
loss function (dashed red) are shown. From the latter 
one can see that only one resonance occurs and thus the 
additionally found solution is physically irrelevant. 



FIG. 6: (Color online) (a) Real part of the dielectric function 
(solid black) and imaginary part of the energy loss function 
(dashed red) of GaAs. (b) The same for InAs. Parameters: 
q = 0.3, <p q = 0°, E z = 0, and n = 5 ■ 10 14 m" 2 . 



B. InAs 

So far we have seen that plasmons in GaAs-based two- 
dimensional hole systems show a pronounced anisotropy 
as manifested, e.g., in the plasmon spectrum and 
damping rate. Unfortunately, as the Rashba contribution 
in GaAs is virtually negligible, see Table [TJ there 
is no direct control of this feature, e.g., by turning 
on/off the anisotropy dynamically by varying an electric 
field. This control, however, is necessary for possible 
applications, e.g., as a plasmon filter or a plasmon field- 
effect transistor, as already proposed for 2DEGs!^ On 
the other hand, as we know that in InAs the Rashba 
parameter and its influence on the energy spectrum is 
large, we expect SOC to give the possibility to modulate 
the plasmon spectrum as well. In principle, one might 
also use materials with an even larger Rashba coupling 
constant than that of InAs, e.g., InSb. However, the 
reason we focus on the former is twofold. First of all, in 
experiments InAs is more popular than InSb and, second, 
the four-band Luttinger Hamiltonian used in this work is 
not sufficient to describe InSb where the more advanced 
8x8 Kane model including electronic conduction bands 
and split off bands needs to be usedP^l 

In Fig. [7] the plasmon energy and the SPC of InAs 
without an electric field is shown for cj)„ — 0° (solid 
black) and <\> q = 45° (solid red). As before, the carrier 
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FIG. 7: (Color online) The solid lines show the plasmon 
dispersion of InAs for an angle orientation of 4> q — 0° (black) 
and (j) q — 45° (red). The dashed lines display the boundaries 
of the single-particle continuum for <j> q = 0° (black) and 
4> q = 45° (red). The electric field is set to E z = and the 
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FIG. 9: (Color online) Electric field dependence of 
the plasmon energy uj q (solid lines) and lifetime being 
proportional to the imaginary part of the dielectric function 
Im {e(q, u q )} (dashed line) for fixed wave vector q — 0.3 and 
two different angle orientations, <f) q = 0° (black) and <f) q = 45° 
(red). The hole density is again set to n = 5 ■ 10 14 m~ 2 
(E F = 0.3 • E ). 



concentration is n = 5 • 10 14 m~ 2 , corresponding to 
Ep = 0.3 • Eq. Region I marks again transitions with 
initial and final states in E^i, while in region II the final 
state is Ehs. Not shown in the plot is the interband 
continuum with final states in En and Ei 2 , respectively. 
Unlike GaAs, the interband and intraband parts merge 
at q 0.6. The anisotropy in the plasmon spectrum 
is less pronounced than in GaAs. The SPC, on the 
other hand, is only weakly direction dependent where the 
main anisotropy arises in the intraband part. Contrary 
to GaAs, the two modes enter the continuum at almost 
equally large momenta q ~ 0.45. 

Up to now, spin-splitting of the energy bands is 
only caused by the Dresselhaus term. The effect of 
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FIG. 8: (Color online) Energy loss function of InAs for 
fixed q = 0.3 and various electric fields E z — (black), 
1.8- 10 7 V/m (green) and 2.8- 10 7 V/m (red) for two different 
angle orientations cj> q — 0° (a) and 45° (b). The hole density 
is set to n = 5 ■ 10 14 m~ 2 (E F = 0.3 ■ E ). 



an increasing electric field, pointing along the growth 
direction, which leads to a finite Rashba contribution in 
Eq. Q , is demonstrated in Fig. [8] for fixed momentum 
q = 0.3 and two different directions <j) q = 0° and <j) q = 
45° . The 5-\ike peak in the case without applied external 
electric field becomes broadened for large enough E z , 
because both the intra- and interband continua expand 
due to the spin-splitting. If the interband SPC is 
sufficiently broad, it finally contains the point of the 
resonant energy. The detailed field dependence of the 
plasmon energy and lifetime is shown in Fig. [9j The 
energy (solid line) is blueshifted once the electric field is 
enlarged whereat the difference between the E z = and 
E z = 3.0 • 10 7 V/m result is about 15 percent. 

More important than the position of the peak is its 
width, which in turn is related to the imaginary part of 
the dielectric function (dashed lines). For <\> q = 0° the 
mode remains undamped if E z < 1.0 • 10 7 V/m. For 
larger fields, the quasiparticles acquire a finite lifetime as 
Im{e(g,w 9 )} is nonzero. Qualitatively the same features 
can be seen for <p q = 45° . However, the important point is 
that the critical field at which the mode enters the SPC, 
E z sa 1.9 • 10 7 V/m, is now much larger. This in turn 
yields two possible applications. First, in a plasmon field- 
effect transistor, excitations in the source can be detected 
in the drain, depending on their lifetime and thus on the 
strength of the applied electric fieldP' Second, changing 
the electric field in the range of 1.0 • 10 7 V/m < E z < 
1.9- 10 7 V/m allows for plasmon filtering since particular 
directions are damped while others are notP 

Finally, we want to remark that as in GaAs, additional 
solutions of the approximate solution Eq. p| can be 
found in InAs. However, as can be seen in Fig.TfJlb), only 
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one of these solutions corresponds to a genuine plasmonic 
mode with a characteristic resonance in the energy loss 
function. 



V. CONCLUSIONS AND OUTLOOK 

We have investigated the dynamical dielectric 
functions of a GaAs and InAs quantum well, 
including spin-orbit interactions of the lowest order 
Dresselhaus and Rashba type, in the four-band Luttinger 
Hamiltonian for the two lowest subbands. For GaAs a 
pronounced anisotropy in the plasmon spectrum and in 
the single-particle continuum occurs. Depending on the 
direction of the incident beam, plasmons are damped or 
long-lived. This opens, in principle, the possibility of 
filtering plasmons with distinct orientations!^ However, 
the main problem in GaAs is the disability of controlling 
this feature due to the negligible Rashba contribution. 

Here InAs seems to be more promising. In fact, it 
turns out that while long-wavelength plasmons in InAs 



do not decay provided the Rashba contribution is small 
enough, the lifetime of the plasmons can be modulated by 
changing the electric field. Furthermore, we have shown 
that the critical field at which a finite damping occurs is 
direction dependent. These features might be useful for 
applications such as a plasmon field-effect transistor^ or 
a plasmon filterPEl 

Finally let us point out that while this study is based 
on realistic standard band parameters,^ the strength 
of the spin-orbit coupling contributions can even be 
enlarged beyond that values by applying strain to the 
sample or using more advanced heterostructures which 
might lead to an enhancement of the effects described in 
this work. 
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